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摘要 : 银河 系 的 邻居 仙女 座 星系 (M31) 是 本 星系 群 中 恒星 质量 最 大 的 成 员 。 观 测 表明 ，M31 正 
与 银河 系 相互 接近 ， 且 两 者 极 有 可 能 在 未 来 发 生 并 合 。 N 体 模拟 作为 研究 星系 并 合 的 有 力 工具 ， 
在 过 去 半 个 世纪 内 取得 了 长 足 的 发 展 ， 也 在 有 关 银 河 系 与 M31 并 合 的 定量 研究 中 发 挥 了 巨大 
的 作用 。 随 着 星系 结构 理论 的 丰富 、 观 测 水 平 的 提高 以 及 N 体 模 拟 软件 技术 的 成 熟 ， 银 河 系 与 
M31 并 合 的 图 像 也 越 来 越 清晰 与 完善 。 简 要 介绍 了 银河 系 与 仙女 座 星系 并 合 模拟 研究 的 工作 流 
程 ， 回 顾 了 各 步 又 所 用 到 的 理论 模型 或 软件 条 件 于 近 几 十 年 内 取得 的 发 展 ， 选 取 各 阶段 代表 性 模 
拟 结果 进行 比较 ， 最 后 给 出 了 未 来 的 研究 趋势 与 改进 方向 。 
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1 引 言 


本 星系 群 (Local Group, LG) 包含 银河 系 (Milky Way, MW)、 仙 女 座 星系 (Messier 31, 
M31)、 三 角 座 星系 (Messier 33, M33)、 麦 哲 伦 云 等 近 百 个 大 小 不 一 、 形 态 各 异 的 星系 ”， 
有 关 本 星系 群 的 研究 对 于 理解 星系 间 相 互 作用 、 星 系 演化 进程 、 宇 宙 学 等 具有 重大 意义 。 其 
中 ，M31 是 本 星系 群 内 恒星 质量 和 光度 最 大 的 星系 ”， 与 银河 系 一 起 构成 了 一 个 可 用 于 研究 
星系 作用 与 结构 演化 的 “天 然 实 验 室 ”。 按 照 宇宙 学 的 等 级 成 团 理 论 ， 星 系 是 并 合 增 长 的 。 
本 星系 群 中 不 仅 有 众多 卫星 星系 的 并 合 ， 未 来 还 可 能 发 生 银河 系 与 M31 这 两 个 大 质量 星系 
的 并 合 。 
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4 期 


主流 观点 认为 中， 本 星系 群 已 经 与 字 


河 系 与 M31 FFE RA AW 


合 的 可 能 性 。 


Gi 


BANK ANE, AF 
响 ， 或 者 说 ， 本 星系 群 内 部 的 星系 已 经 构成 一 个 引力 束缚 系统 。 
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HS BINS Ziyi (Hubble flow) 的 影 
因此 在 万 有 引力 作用 下 ， 银 


观测 结果 显示 ，M31 确实 具有 朝向 和 


河 系 的 


径 向 速度 ，Binney 和 Tremaine” 于 1987 年 给 出 该 相对 速度 约 为 -120lm . s-1; 而 近 些 年 
给 出 的 最 新 径 向 速度 略 低 ， 约 为 -109km.s-1”。 由 于 切 向 速度 对 两 者 的 相对 运动 也 存在 较 


影响 ， 


年 间 ， 对 M31 
合 事 件 的 描述 也 愈加 准确 。 


化 、 星 系 并 合 的 重要 工具 。 
合 的 可 能 性 后 ， 能 较 好 


关于 银河 系 与 M31 并 合 的 定量 结论 ， 如 两 者 完成 并 合 的 时 间 、 


也 在 不 断 更 新 优化 ， 这 种 更 新 3 


位 置 等 。 银 河 系 与 M31 并 合 的 过 程 图 像 


(1) 星系 结构 理 


可 以 用 更 合适 的 数学 模型 来 描述 星系 的 各 个 组 成 成 分 (暗物质 坚 、 
测 水 平 的 提高 ， 星 系 总 质量 (这 里 一 般 指 位 力 质 量 )、 星 系 各 


长 、 标 高 、 半 光 半 径 等 ) 等 参数 愈加 明晰 


在 人 们 根据 观测 数据 以 及 理论 分 析 推 断 出 外 
苗 述 星系 具体 结构 与 星体 运动 的 入 体 模拟 工作 也 逐步 开展 ， 给 出 了 
并 合 完成 后 太阳 可 能 所 处 的 
要 基于 三 个 方面 ; 
论 的 丰富 ， 人 们 在 构建 星系 模型 时 能 够 选择 更 多 种 类 的 分 布 轮廓 来 拟 合 ” ， 


ad (如 机 
BET ike A 


HR AR AY 


因此 在 M31 的 自行 被 测定 之 前 ， 仅 仪 通过 径 向 速度 无 法 完全 确定 它 是 否 会 与 外 
系 并 合 ， 不 过 人 们 可 以 根据 本 星系 群 提供 的 其 他 
星系 的 轨道 等 ) 间接 推断 出 M31 切 向 速度 的 可 能 区 间 ，3 


自行 的 测量 已 经 取得 初步 的 结果 ”， 相 关 工 作对 银河 系 与 仙女 座 星 系 未 来 并 


BRE 


河 


寺 构 特征 、 某 些 矮 
的 模拟 工作 。 过 去 10 


自 20 世纪 70 年 代 Toomre A 和 Toomre J" 利用 计算 机 粒子 数值 模拟 研究 星系 间 相互 
以 来 ， 数 值 模拟 技术 (包括 N 体 动 力学 模拟 和 流体 模拟 ) 不 断 发 展 } 


成 为 研究 星系 演 


河 系 与 M31 存在 并 


数值 模拟 中 特殊 问题 的 经 验 积累 。 


值得 指出 的 是 ， 所 谓 “ 理 想 星系 并 合 ” 的 模拟 通常 注重 研究 不 同 物 天 


角 动 量 等 ) 对 并 合 结果 的 影响 ， 并 聚焦 于 3 


他 并 合 系统 所 不 具备 的 观测 限制 。 


本 文 将 回顾 几 十 年 来 银河 系 与 M31 并 合 模拟 的 下 
系 并 合 N 体 模拟 的 操作 流程 展开 ， 简 要 介绍 
旨 在 让 读者 能 够 利用 文章 中 介绍 的 知识 与 工具 自主 实现 4 


发 展 ， 
因此 本 文 也 可 以 作为 星系 并 合 N 体 模 拟 了 


ERA ARE); (2) 观 
分 质量 占 比 、 星 系 尺 度 ( 标 


m= M31 的 运动 学 参数 愈加 精确 "一 ;， (3) N 体 
模拟 软件 技术 的 成 熟 旦 号， 不 仅 包 括 程序 的 完备 化 与 便捷 化 ， 也 包括 算法 的 优化 以 及 处 理 


步骤 中 所 


究 进 展 。 值 得 强调 的 是 ， 文 章 按照 星 


[ 作 的 参考 手册 。 


Poke 


5 


发 展 与 观测 数据 更 新 对 上 述 选择 的 影响 ， 第 3 章 将 介绍 模拟 所 需 初始 条 伯 


要 包括 单个 星系 初始 条 件 生 成 器 的 开发 ， 


参数 (如 质量 比 、 
合 对 星系 中 心 超大 质量 黑洞 活动 性 的 影响 ， 这 方 
面 工作 的 综述 可 见 参考 文献 Im。 而 银河 系 与 M31 的 并 合 ， 作 为 星系 并 合 的 一 个 特例 ， 有 其 


=I 


到 理论 模型 与 软件 条 件 的 更 新 与 
MAS M31 的 并 合 模 拟 ， 


2 章 将 介绍 银河 系 与 M31 基本 参数 的 选取 以 及 物质 分 布 轮廓 的 选择 ， 表 明基 础 理论 
文件 的 生成 ， 


ry 


以 及 人 们 对 星系 空间 关系 和 速度 关系 的 认 知 进展 ; 
第 4 章 将 介绍 N 体 模 拟 所 用 软件 ， 简 述 相关 软件 的 更 新 以 及 问题 处 理 经 验 的 积累 ， 第 5 章 
将 展示 所 选 代表 性 工作 的 模拟 结果 ; 第 6 章 将 进行 总 结 并 提出 有 价值 的 未 来 研究 方向 。 
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2 星系 模型 的 构建 


星系 包含 暗物质 、 恒 星 、 


气体 、 中 心 黑洞 


tH A 


星 ， 其 中 恒星 主要 存在 于 星系 盘 与 核 球 之 


41 卷 


N 体 模拟 重点 考虑 的 是 暗物质 党 与 
Ph。 为 了 进行 银河 系 与 M31 的 数值 模拟 工作 ， 


=} 


WR 


球体 内 的 平均 密度 为 200 倍 宇宙 临界 密度 )、 自 
BAAS). RRB (He. ARE. 


T RE 
(半径 、 厚 度 等 )。 


首先 要 确定 星系 的 基本 参数 ， 主 要 包含 : 星系 质量 、 
程度 ， 定义 为 c= r200/Ts> 其 中 C NEREZ, Ts 为 标 长 暗物质 晕 在 半径 为 7200 
eB (spin parameter, 


HES) Joe E BY) ae: RS EM ( 


述 星系 中 物 


H 


核 球 等 ) 质量 、 


一 般 用 字母 入 表 
星系 各 组 分 尺度 


不 过 ， 以 上 基本 参数 尚 不 足以 勾勒 出 星系 的 形态 结构 ， 还 需要 知道 质量 在 星系 中 是 如 何 
分 布 的， 这 就 需要 选取 合适 的 分 布 函 数 对 各 个 成 分 进行 拟 合 。 


值得 一 提 的 是 ， 为 了 使 所 构建 的 星系 更 加 真实 可 信 ， 


某 些 参数 细节 ， 试 图 使 其 自转 曲线 


也 为 模拟 得 弄 


a 


与 观测 到 的 星系 自转 曲线 在 某 一 3 


总 质量 固定 上 


体 结论 (如 银河 系 与 M31 


| 的 某 些 具 
本 章 将 回顾 近 些 年 来 人 们 对 和 


的 银河 系 与 M31 模型 。 
2.1 FAAS M31 的 质量 


前 文 提 到 的 星系 参数 有 若干 个 ， 而 此 处 着 重 讨论 质量 ， 一 方面 是 因 
质量 显然 是 所 有 参数 中 对 两 者 未 来 演化 影响 最 大 的 ， 另 一 方面 ，] 


确定 星系 质量 后 调整 其 他 参数 (如 


接着 介绍 几 种 星系 数值 模拟 中 经 典 的 分 布 轮廓 ， 最 后 简 自 


Bx Ae j 
REE 


百 星 质量 占 总 质量 的 


近 真 实情 况 ， 所 以 这 些 参数 并 不 
选取 的 。 下 面 我 们 回顾 人 们 对 
数值 。 


JE re FA WL 


群 总 质量 ， 第 二 个 环节 是 确定 和 


其 核心 思想 是 


认为 4 


膨胀 而 相互 远离 ， 随 后 这 种 远离 的 运动 又 发 4 


河 系 与 M31 的 质量 比 。 对 本 与 
依靠 计时 推理 法 (timing argument)。 这 一 技术 最 早 于 1959 年 由 
河 系 与 M31 在 宇宙 致密 早期 阶段 诞生 时 非常 接近 ， 接 着 由 于 宇宙 
E 反 转 ， 两 者 在 引力 作用 下 相互 靠近 ， 如 


1 得 到 ， 而 是 根据 


期 人 们 对 于 银河 系 与 M31 质量 的 估计 分 为 两 个 环节 ， 第 一 个 环节 是 确定 本 


E 如 前 文 所 述 ， 人 
比例 等 ) 使 其 
些 分 析 计 算 与 试探 | 
Li] RAI M31 质量 的 认 知 演化 以 及 在 相关 模拟 工作 中 选取 的 


的 情况 下 ， 人 们 会 调整 
# 径 范围 内 尽量 符合 ， 这 


并 合 后 太阳 将 处 于 何 处 ) 增添 说 服 力 。 
河 系 与 M31 质量 的 选取 随 观测 数据 更 新 而 产生 的 变化 ， 
展示 相关 数值 模拟 工作 中 构建 出 


为 两 个 星系 各 自 的 
门 会 在 
旋转 曲线 接 
生 模 拟 


星系 


Lo hy 


EE 系 群 总 质量 的 估计 主要 
Kahn 和 Waltjer” JF fil, 


今 它 们 的 轨道 已 经 完全 符合 开 普 勒 定律 。Binney 和 Tremaine” F 1987 年 的 工作 中 应 用 


了 该 技术 ， 在 确定 银河 系 与 M31 We 
量 应 大 于 3 x 10! Mo; 同时 他 们 也 估算 出 和 


群 的 总 质 


将 于 4Ga 以 后 到 来 。 当 然 这 只 是 一 个 粗略 的 结果 。 
Ri AD M31 这 两 个 星系 贡献 ， 但 也 
的 贡献 不 可 


为 其 主要 由 外 
medium, IGM) 对 总 质 


= 
里 


1 更 一 般 地 说 是 暗物质 晕 的 质量 分 布 轮廓 。 坚 聚 身 


E 离 、 径 向 相对 速度 、 宇 宙 年 龄 后 
河 系 与 M31 的 下 一 次 近 心 点 时 刻 
系 群 总 质量 后 ， 人 们 认 
理论 认为 星系 


A 


在 给 出 本 星 


， 分 析出 本 星系 


际 介 质 (intergalactic 


忽略 ， 这 里 为 简单 说 明 原 理 而 暂 不 考虑 之 。 利 用 


度 只 是 其 中 较 常见 的 一 种 参数 ， 常 用 于 本 文 2.2.2 节 提 到 的 NFW 轮廓 。 


4 期 


Ruw = (4.5 加 1.0) kpc (1 
1991 年 给 出 )、 RM31 = (5.8 lee 


周 宇 凡 ， 等 : 


计算 出 M31 与 
Tully-Fisher 关系 


河 系 质量 比 : 


ea . 
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星系 旋转 曲线 平坦 部 分 的 数值 以 及 B 波段 观测 到 的 星系 盘 标 长 : Vw = 220km-s7!, 
Freeman” 在 1987 EAH) Vegi = 260 km - s7! (由 Braun” # 


0.3) kpc (由 Walterbos 和 Kennicutt ™ 在 1988 年 给 出 ) 可 以 
Mm3ı/Mmw = (RusiVi431)/(RuwVaiw) © 1.9。 或 者 根据 


Ma«V® , (1) 


Ma=3~ 4, Af Mmı/Mumuw © 1.5 ~ 2。 上 述 讨 论 详 见 Raychaudhury 和 Lynden- 


Bell ™ 1989 年 的 了 


[ 作 。 


几乎 是 最 早 的 利用 N 体 模 拟 完整 描述 银河 系 与 M31 未 来 并 合 过 程 的 工作 


Dubinski 


a AF 1996 年 发 表 的 工作 (下 文 简称 D96)， 根 据 当 时 由 计时 推理 法 得 到 的 结果 选取 了 两 
个 本 星系 群 总 质量 可 能 值 : Mig jow = 1.6 x 102 Mos MLG nigh = 5-3 x 1012 Mea。 同时 也 根 
据 前 文中 的 有 关 计 算 确 定 银 河 系 与 M31 的 质量 比 Muw : Musi © 1 :2。 因 此 该 工作 中 两 组 


模拟 采用 的 银河 系 与 M31 质量 分 别 为 AMwwl = 5.0 x 104 Mo. Mms1,1 = 1.1 x 10” Mo 以 
及 Mmw.2 =i 1012 Mo. MM31,2 =3.6x 1012 Mo» 
2002 Æ, Klypin 等 人 采用 另 一 种 方法 给 出 了 银河 系 与 M31 的 质量 。 他 们 从 一 些 基 


础 理论 出 发 构建 一 系列 模型 ， 这 些 模型 的 总 质量 、 质 量 分 布 、 旋 转 曲 线 、 具 体 结构 形态 等 


各 不 相同 ;接着 根据 星系 的 观测 数据 给 定 一 系列 限制 条 件 (比如 距 星 系 中 心 一 定 范围 内 所 包 


含 的 质量 为 多 少 )， 调 整 各 模型 
束 条 件 兼容 ， 便 认为 该 模型 是 合理 的 


的 参数 细节 以 进行 拟 合 ， 经 反复 试验 ， 当 某 个 模型 成 功 与 约 


该 工作 并 不 意欲 寻找 “最 好 ”的 模型 与 参数 ,所 以 最 


终 给 出 了 一 系列 符合 所 选取 约束 条 件 的 银河 系 与 M31 模型 ， 而 作者 认为 其 中 最 “合理 ”( 指 


星系 盘 质 量 、 曙 聚集 度 等 参数 在 各 自 可 能 的 取 值 区 间 内 最 为 适中 ) 的 两 个 星系 模型 质量 分 别 


为 ， Mrw = 1.0 x 10!2 Mo, Mm3ı = 1.6 x 101? Ms。 该 数值 在 Cox 和 Loeb™ 2008 年 的 NN 
体 模拟 工作 (下 文 简称 C08) 以 及 Schiavi 等 人 ”2020 年 的 N 体 模 拟 工作 (下 文 简称 S20) 中 


被 采用 。 


2012 年 ，van der Marel 等 人 根据 哈 勃 空间 望远镜 (Hubble Space Telescope, HST) 对 
运用 计时 推理 法 ， 结 合 前 人 文献 中 的 一 些 讨 论 ， 得 出 新 的 本 星系 群 总 质 


M31 的 观测 数据 ， 
量 : Mic = (3.17 


工 0.57) x 10! 


2 Moo 21 世纪 以 来 ， 在 van der Mare 这 篇 文章 发 表 前 ， 也 


存在 其 他 利用 计时 推理 法 估计 本 星系 群 质量 的 工作 ( 详 见 参考 文献 回 中 的 介绍 及 其 参考 文 
献 )， 它 们 给 出 的 质量 更 大 ， 大 于 人 们 普遍 认为 的 银河 系 与 M31 质量 之 和 ， 因 此 N 体 模 拟 


工作 中 必须 添加 星系 际 介 质 以 


使 模拟 包含 的 总 质量 与 本 星系 群 总 质量 相符 ， 而 该 较 小 的 本 


ERRI. EET, 


星系 群 质 量 估 计 结 果 出 现 后 ,多 


相 比 ， 其 认为 最 新 的 结果 接近 了 
] Muw = Musi = 1.5 x 10 Mo 的 星系 质量 进行 了 N 体 模拟 工作 (下 文 简称 V12)。 
随 着 越 来 越 多 观测 数据 (如 Gaia 数据 ) 的 发 布 ， 银 河 系 与 M31 质量 的 研究 结果 也 越 


这 一 工作 对 了 


河 系 与 M31 的 并 合 模拟 工作 中 星系 际 介质 便 不 再 是 必须 考 
F 银河系 与 M31 的 质量 比 也 给 出 了 不 同 的 结果 ， 与 结果 1 :2 


1 : 1。 于 是 van der Mared 等 人 四 在 同系 列 另 一 篇 文章 中 采 


来 越 丰富 。 如 2019 Æ Watkins 等 人 "根据 Gaia DR2 数据 ， 结 合 其 他 工作 分 析出 银河 系 
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质量 应 为 Muw = 1.54437) x 102 Mo, X5 V12 中 所 采用 的 银河 系 质量 相符 。 然 而 银 
河 系 与 M31 的 质量 比 并 没有 一 个 公认 的 结果 ，1 : 2 和 1 :1 RRA, Me 
能 性 。Patel 等 人 ”于 2017 年 的 一 项 研究 大 麦哲伦 云 与 M33 轨道 历史 的 工作 中 ， 对 银河 
R5 M31 各 构建 了 两 个 模型 ， 其 中 Mumwi = 1.0 x 102 Mo. Mw. = 1.5 x 10! Mos 
Mm311 = 1.5 x 101? Mo, Mm31,2 = 2.0 x 10!2 Mo, van der Marel 等 人 "于 2019 年 再 次 进 
行 的 N 体 模拟 工作 (下 文 简称 V19) 中 采用 上 述 1:2 和 1 :1 两 种 质量 比 模 型 分 别 进行 了 
模拟 。 

对 于 银河 系 与 M31 质量 的 研究 与 争论 也 从 未 终结 。2022 年 底 ，Patel 和 Mandel” iE 
据 4 个 卫星 星系 (M33. NGC185. NGC147,. IC10) 的 轨道 数据 ， 给 出 M31 的 质量 可 能 大 
Mms: = 2.817978 x 1012 Mo 甚至 Musi = 2.995522 x 1012 Mo， 这 比 之 前 的 估计 值 要 高 很 
多 。 这 也 为 未 来 的 模拟 工作 提供 新 的 可 能 参数 。 

2.2 ”物质 分 布 轮廓 

HARI. BAG (恒星 盘 与 气体 盘 )、 核 球 等 成 分 组 成 ， 各 成 分 所 包含 的 物质 分 
布 形式 一 直 是 人 们 感 兴趣 的 话题 。 星 系数 值 模拟 工作 中 ， 由 于 要 在 计算 机 上 构建 星系 模型 ， 
合理 的 物质 分 布 函数 表达 式 便 显 得 尤为 重要 。 过 去 几 十 年 间 ， 经 过 大 量 工作 的 试验 与 讨论 ， 
不 少 物 质 分 布 轮廓 被 提出 。 有 具体 的 研究 过 程 与 发 展 历程 十 分 复杂 ， 这 里 不 加 歼 述 ， 仅 介绍 N 
体 模 拟 中 最 常用 的 几 种 分 布 。 

2.2.1 Hernquist 轮廓 (适用 于 核 球 与 暗物质 晕 ) 

人 们 很 早 就 认识 到 ， 椭 圆 星系 以 及 星系 核 球 的 面 亮度 分 布 可 以 被 de Vaucouleurs ™ 的 

RL 定律 很 好 地 描述 : 
R 1/4 
(z) - | 


其 中 ， 玉 为 距 星 系 中 心 的 投影 半径 ，R。 为 有 效 半径 或 者 半 光 半径 ，J 为 表面 亮度 。 尽 管 该 
函数 形式 能 很 好 地 描述 椭圆 星系 以 及 星系 核 球 ， 它 所 对 应 的 空间 密度 分 布 函数 却 没有 解析 
解 下 ， 也 就 难以 用 于 理论 分 析 ， 所 以 寻找 一 个 便于 解析 且 投 影 面 亮 度 轮 廓 类 似 于 R14 的 空 
间 密 度 函 数 就 显得 尤为 重要 。 经 过 一 些 尝试 ， 最 终 Hernquist 轮廓 ”问世 ， 其 密度 轮廓 为 : 


lg Foal = 一 3.331 


_Ma 1 
— nr (r+a)3 


p(r) (3) 


其 中 ，p 是 密度 ，MM 是 质量 ，a 是 标 长 ，” 是 距 星系 中 心 的 三 维 半径 。 这 一 轮廓 在 星系 数值 
模拟 中 被 广泛 使 用 ， 通 常 被 用 来 描述 核 球 以 及 暗物质 坚 的 物质 分 布 。 
2.2.2 NFW 轮廓 (适用 于 暗物质 军 ) 

Navarro 等 人 "提出 的 NFW (Navarro-Frenk-White) 密度 轮 廊 表达 式 为 : 


pr) _ bc 
Perit (r/rs\1+r/rs)? : (4) 
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其 中 ，p 是 密度 ，6。 是 无 量 纲 化 的 特征 密度 ，rs 是 标 长 ，” 是 距 星 系 中 心 的 三 维 半径 ，perit 
是 闭 全 而 界 密度 a 
Perit = er ) (5) 


KH, H 是 哈 勃 常数 ，G 是 万 有 引力 常数 。 

NEW 轮廓 一 般 仅 被 用 于 描述 暗物质 坚 的 物质 分 布 ， 其 与 Hernquist 轮廓 最 大 的 不 同 
在 于 极 远 处 的 渐进 行为 : 当 7 > r 时 ，NFW 轮 廊 函数 趋 于 -3 而 Hernquist 轮廓 函数 趋 
于 4。 所 以 ， 对 于 同等 质量 的 星系 暗物质 晕 ， 采 用 NFW 轮廓 描述 会 使 其 分 布 得 比 采 用 
Hernquist HERE “mP”, Haapse Ke, NEW 轮廓 在 应 用 时 通常 要 搭配 一 个 人 
为 设置 的 截止 半径 。 

如 果 模 拟 只 对 宏观 的 结果 (如 两 星系 最 终 并 合 的 时 间 ) 感 兴趣 而 不 在 乎 星系 内 部 细节 ， 

么 星系 盘 与 核 球 选择 何 种 轮廓 影响 并 不 大 ;而 暗物质 晕 占 据 着 星系 高 于 90% 的 绝 大 部 分 

质量 ， 决 定 着 星系 延伸 至 何 处 ， 所 以 暗物质 晕 的 分 布 轮廓 对 模拟 结果 影响 明显 。Springel 等 
ATHE: Hernquist 轮廓 在 靠近 中 心 的 部 分 与 NEW 轮廓 极其 相似 ， 而 因为 其 轮廓 曲线 在 
外 部 下 降 得 比 NEW 更 快 所 以 总 质量 收敛 ， 于 是 采用 Hernquist 轮廓 描述 暗物质 晕 时 便 不 需 
要 人 为 设置 截止 半径 ， 另 外 ，Hernquist 轮 廊 具有 可 解析 的 分 布 函数 形式 ， 这 样 在 分 析 时 比 
NEW 轮廓 方便 得 多 。 综 上 所 述 ，Hernquist 轮廓 比 NFW 轮廓 更 适合 在 N 体 模拟 中 描述 暗 
物质 坚 
2.2.3 ”指数 轮廓 (适用 于 星系 盘 ) 

指数 轮廓 表达 式 为 : 


Y 
TL 


o(R) = expl -R/h) ， (6) 
Kt, p NE, MAME, h 为 标 长 ，R 是 距 星 系 中 心 的 二 维 半 径 ， 这 一 轮廓 经 常用 来 
描述 星系 盘 的 物质 分 布 。 由 于 星系 盘 并 不 是 一 个 球 对 称 系统 ， 上 述 表达 式 仅 能 描绘 径 向 分 
布 ， 如 果 要 完整 表示 三 维 空间 中 盘 的 物质 分 布 ， 还 需要 加 上 一 个 纵向 成 分 ， 比 如 sech? 项 ， 
则 指数 轮廓 变 为 : 


M of Z 
p(R, z) = mar S (=) exp(—R/h) ， (7) 


其 中 ，z 为 距离 盘面 的 高 度 ，zo 为 标高 。 
2.2.4 Miyamoto-Nagai Æ (适用 于 星系 盘 ) 

为 了 描述 轴 对 称 且 有 厚度 的 星系 盘 ，Miyamoto-Nagai 得 "被 引入 ， 这 是 构建 星系 模型 
时 常见 的 一 种 分 布 ， 它 的 势 函 数 为 : 


&(R,z) = 


GM 
yR? + (Rat (e+e? 
其 中 ，5 为 势 ，G 为 万 有 引力 常数 ，M 为 质量 ，R 为 距 星 系 中 心 的 二 维 半径 ，z 为 距离 盘 
面 的 高 度 ， Ra 5 za 为 标 长 与 标高 。 

随 着 Ra 与 za 取 值 的 变化 ， 该 势 函数 可 以 描述 从 薄 盘 到 球 等 多 种 模型 ， 模 型 无 奇 点 且 
函数 简洁 ， 模 拟 时 的 计算 成 本 比 指数 盘 等 复杂 的 轮廓 更 低 。 


(8) 
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2.2.5 Plummer 球 (适用 于 核 球 ) 
Plummer 球台 的 势 函 数 为 : 


GM 


J (9) 


P(r) = 一 


其 中 ，5 为 势 ，G 为 万 有 引力 常数 ，M 为 质量 ,7 为 距 星系 中 心 的 三 维 半径 ，a 为 标 长 。 该 
模型 一 般 用 于 拟 合 星系 核 球 ， 在 某 些 特别 的 工作 中 对 于 不 需要 分 辨 具体 细节 的 星系 也 可 以 
使 用 该 模型 。 比 如 研究 卫星 星系 绕 大 星系 运行 的 轨道 时 ， 不 需要 分 辨 卫星 星系 内 部 诸如 星系 
得 、 核 球 之 类 的 具体 细节 ， 这 时 就 可 以 用 一 个 Plummer 球 来 表示 整个 卫星 星系 。 
2.3 ”银河 系 与 M31 模型 示例 

在 确定 了 星系 参数 与 各 成 分 分 布 轮廓 后 ， 星 系 模型 也 就 基本 构建 完成 。 下 面 以 C08 为 
例 介绍 其 中 选用 的 银河 系 与 M31 模型 。 此 模型 中 的 星系 参数 由 前 文 提 到 的 Klypin 等 人 号 在 
2002 年 的 工作 给 出 ， 其 也 在 S20 中 被 采用 。 

表 四 展示 了 C08 与 S20 选取 的 星系 参数 ， 其 中 暗物质 晕 和 核 球 采用 Hernquist 轮廓 ， 


星系 盘 采 用 指数 轮廓 。 


#1 Cos 与 $20 中 采用 的 星系 模型 


参数 MW M31 
V200 / km:s-! 145 170 
M200 / 10° Mo 1.0 1.6 
c 12 12 
入 0.031 0.036 
ma 0.041 0.044 
Ra / kpc 2.2 3.6 
f 0.3 0.3 
Mp 0.008 0.012 
a / kpe 0.4 0.7 


YE: 选 自 C08 中 的 表 1。 第 2 一 10 行 分 别 指 位 力 速度 、 位 力 质量 、 晕 聚集 度 、 自 转 参数 、 星 系 盘 质量 
分 数 (星系 盘 质 量 占星 系 总 质量 的 比例 )、 星 系 盘 标 长 、 星 系 盘 中 气体 比例 (严格 意义 上 N 体 模拟 不 包 
含 气 体 ， 此 处 也 用 到 一 部 分 流体 模拟 )、 核 球 质量 分 数 ( 核 球 质量 占星 系 总 质量 的 比例 )、 核 球 标 长 。 


不 过 ， 星 系 模型 并 不 一 定 要 完全 按照 这 一 示例 ， 星 系 参数 并 不 一 定 是 表 中 所 列 ， 读 者 可 
以 查阅 前 文 提 到 的 另外 三 个 工作 进行 对 比 。 另 外 ， 物 质 分 布 轮 亡 也 可 以 根据 需要 调 
整 ， 比 如 V12 中 的 核 球 即 采 用 R14 HERR, Patel 等 人 对 银河 系 与 M31 RMD ER 
NFW 轮廓 ， 对 星系 盘 采 用 Miyamoto-Nagai 盘 。 但 近 些 年 的 N 体 模拟 中 ,“Hernquist 上 暗 
物质 晕 + 指数 盘 十 Hernquist 核 球 ” 似 乎 是 一 个 比较 “优质 ”的 选择 ， 本 文 2.2.2 节 与 3.1 节 


的 相关 内 容 对 此 给 出 解释 。 
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3 星系 初始 条 件 文件 的 生成 


银河 系 与 M31 的 到 
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LE 论 模 型 构建 成 


Iya, fi 


究 人 员 需 要 在 计算 机 中 将 它们 创建 出 来 ， 这 


是 一 个 比较 耗费 精力 的 步 又。 随 着 模拟 工作 “程式 化 ”程度 越 来 越 高 ， 人 们 迫切 希望 能 


一 个 工具 


只 要 输入 参数 就 可 以 生成 一 个 由 大 量 粒 子 表示 的 星系 。 近 20 年 来 (尤其 是 最 


近 10 年 )， 几 款 初始 条 件 (initial condition, IC) 文件 生成 器 公开 发 布 ， 大 大 简化 了 该 步骤 的 


流程 。 


然而 即使 如 此 ， 这 一 步骤 也 并 没有 那么 简单 。 NN 体 模拟 所 
[ 作 是 极 
平衡 稳定 (仅仅 将 一 
因此 每 个 粒子 的 初始 位 置 和 初始 速度 都 应 当 恰 到 好 处 ) 的 星系 ， 
傅 形 下 ， 生 成 初始 条 作 


H 


起 着 决定 性 作用 


这 个 星系 必然 是 不 稳定 
IC 文件 生成 器 可 能 


生成 一 个 好 的 IC SCH 
的 出 现 使 操作 更 加 便捷 ， 但 如 果 要 生成 


的 ， 


万 然 要 运算 很 长 时 间 ， 
之 后 进行 模拟 的 时 间 更 久 。 此 外 ，IC 文件 生成 器 往往 存在 一 定 的 局 限 性 ， 


可 能 是 针对 某 种 特定 的 情境 所 编 
文件 生成 器 的 输出 文件 
还 需要 人 工 将 两 者 合并 成 一 个 新 IC 文件 ， 输 出 文人 


部 分 IC 


IC 文件 后 ， 


成 器 也 只 是 将 合并 这 一 步骤 


, 
= 
J» 


内 置 。 而 在 合 3 


F 对 于 整 


个 了 


N 


重要 的 。 


质量 很 高 、 


=N 


pas 


H 


ERAGE 


至 在 某 


EAH 


JEK } 


a 


重要 ， 随 着 观测 的 进步 ， 这 些 数据 也 在 不 断 更 新 、 完 善 。 


本 章 将 首先 介绍 几 


改变 。 


th EA 
球星 


3.1 ”星系 初始 条 件 生成 器 
2014 年 ，Yurin 与 Springe ”发 布 了 一 款 IC 文件 生成 器 一 一 GALIC (GALaxy Initial 


Conditions). 


的 IC 文人 


F 对 模拟 的 结果 


EH 


然 IC 文件 生成 器 


=I 
里 


F} 并 不 一 定 完全 符合 使 
只 可 包含 一 个 星系 ， 这 样 在 分 别 生 成 银河 系 与 M31 的 
F 包 含 多 个 星系 的 IC 文件 生 
时 ， 银 河 系 与 M31 的 相对 位 置 和 相对 速度 尤为 


RIC 文件 生成 器 ， 分 析 它 们 的 功能 与 缺陷 ， 之 后 介 纤 
体 模 拟 工 作 中 ， 人 们 对 银河 系 与 M31 相对 位 置 和 相对 速度 的 型 


解 随 观测 数据 更 新 而 


粒子 堆砌 在 一 起 ， 


F 所 花 的 时 间 比 
因为 它 在 编写 时 
] 者 的 要 求 。 


相关 NN 
产生 的 


GALIC 是 一 个 从 初始 参数 条 件 出 发 ， 使 用 迭代 方法 构建 N 体 星系 初始 条 件 文件 的 开 


源 代码 。 它 从 初始 粒子 速度 出 发 ， 将 1 
介 值 函数 进 


差异 ， 用 一 个 


无 碰撞 玻 尔 效 曼 方程 的 稳定 解 。 对 于 给 定 的 密度 结构 ， 由 了 


PA ay 


的 约束 来 丰富 价值 函数 ， 进 而 得 到 最 到 


RAX 
RIC 文件 的 强大 


H 


LaF o 


其 改善 了 拟 合 结果 时 


才 会 被 保留 。 


GALIC 创建 出 的 星系 包含 暗物质 曼 、 
物质 晕 和 核 球 采 用 Hernquist 轮廓 ”， 恒 星 盘 采用 指数 轮廓 。 这 一 轮廓 组 合 也 是 Springel 在 


Rs 


一 系列 工作 中 所 采 
被 广泛 接受 和 使 用 
在 下 载 并 编译 


好 GALIC 后 


mj 


, 


ZN 


恒 


EH 


要 在 参数 文人 


FP 输 入 星 


粒子 轨道 产生 的 时 间 平 均 密 度 响应 与 初始 密度 分 布 的 


行 表征 ， 通 过 和 迭代 调整 初始 速度 从 而 最 小 化 该 价值 函数 ， 继 而 得 到 


六 分布 函 数 通 


常 不 唯一 ， 所 以 代码 


想 的 解 。 速 度 调整 则 采 月 
代码 中 使 用 的 方法 收敛 迅速 且 较 为 灵活 ， 是 构建 


星 盘 、 核 球 三 个 成 分 (无 法 4 


JAY, Blt “Hernquist 暗物质 党 十 指数 盘 十 Hernquist 核 球 ” 模 式 的 星 


系 基 本 参数 


随机 过 程 进行 ， 但 


H 
Fe 


成 气体 )， 其 中 暗 


系 


以 及 星系 各 成 分 粒子 
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数 ， 就 可 以 利用 它 创建 一 个 孤立 星系 的 IC 文件 ， 并 且 GALIC 支持 并 行 计算 以 节约 制作 IC 
文件 的 时 间 。 不 过 GALIC 存在 很 大 的 局 限 性 : 首先 ， 物 质 分 布 轮廓 被 限定 ， 使 用 者 无 法 选 


tT 


XH 


择 其 他 轮廓 ， 其 次 ， 生 成 的 IC 文件 无 法 包含 气体 、 中 心 黑 洞 等 成 分 。 需 要 注意 : 虽然 参数 
里 有 设 定 黑洞 质量 的 选项 ， 但 公开 发 行 的 版 本 是 


漏 写 了 生成 黑洞 的 部 分 ， 即 使 填写 了 相 


关 参 数 也 并 不 会 真 的 有 黑洞 生成 。 如 有 可 能 ， 使 用 者 须 自己 对 相关 代码 做 出 修改 。 


同一 年 ， 另 一 款 功能 更 为 强大 的 IC 文件 4 


E 成 器 DICE™ (Disk Initial Conditions 


Environment) 发 布 ， 它 的 自由 度 比 GALIC 高 很 多 。DICE 可 以 任意 添加 星系 成 分 ， 且 有 


16 种 较为 常见 的 分 布 轮廓 ， 因 此 星系 各 成 分 质量 分 布 的 选项 大 大 增加 ， 这 使 其 可 生成 的 星 


系 形态 极为 丰富 (甚至 不 一 定 是 星系 ): WEA MAER FOR. IOLEE Bite SFE 
DICE 内 置 了 合并 两 个 孤立 星系 IC 文件 的 操作 ， 使 用 起 来 更 加 便捷 。 不 过 DICE 也 存在 一 


定 缺 陷 : 由 于 自由 度 太 高 ， 可 供 设置 的 参数 太 多 ， 宛 长 的 参数 文件 
的 IC 文件 中 可 能 额外 存在 研究 者 并 不 感 兴趣 的 参数 (如 恒星 的 年 龄 、 


除了 公开 发 行 的 IC 文件 生成 器 ， 还 存在 一 些 未 公开 的 代码 。 


F 让 人 眼花 综 乱 ， 并 且 生 成 


金属 丰 度 )。 
如 许多 研究 者 使 用 了 


Springel 的 MakeDisk (或 MakeGalaxy) 代码 及 其 修改 版 ， 这 些 代码 使 用 广泛 但 必须 得 到 


Springel 或 者 合作 者 的 六 
收缩 : 重子 冷却 并 聚集 


F 可 。 值 得 注意 的 是 ， 一 些 代码 对 NFW 轮廓 的 暗物质 晕 考 虑 了 绝热 


在 暗物质 坚 9 


一 般 来 说 ， 生 成 孤立 星系 的 IC 文件 后 ， 需 要 让 它 独 自演 化 一 段 时 间 。 一 方面 检验 稳定 
性 ， 研 究 其 结构 是 否 会 随时 间 而 改变 ， 另 一 方面 对 IC 文件 进行 弛 豫 化 处 理 ， 使 其 内 部 分 布 
更 加 自然 。 如 果 生 成 的 IC 文件 质量 足够 高 ， 这 一 步骤 可 能 并 不 必要 ， 对 最 终 模 拟 结果 的 影 
响 不 大 。 但 是 在 茶 些 复杂 的 模拟 工作 中 ， 这 不 可 或 缺 。 这 一 步骤 并 未 被 IC 文件 生成 器 内 置 ， 
是 因为 其 通常 需要 大 量 的 计算 资源 ， 导 致 计算 时 间 显 著 增加 ; 另外 ， 如 果 将 这 一 步 又 直接 置 


FPF 心 ， 这 一 过 程 将 会 加 深 上 暗物质 党 中 心 区 域 的 重力 势 阱 。 
绝热 收缩 效应 对 星系 的 轮廓 尤其 是 速度 轮廓 影响 显著 ， 经 其 修正 后 的 NFW He E ES 
Ea 


A IC 文件 生成 器 ， 那 么 每 次 需要 重新 生成 IC 文件 时 都 要 执行 完整 演化 过 程 ， 使 程序 变 得 


繁琐 且 难 以 维护 。 


3.2 M31 的 运动 学 参数 
分 别 生 成 银河 系 与 M31 的 IC 文件 后 ， 就 要 根据 两 者 的 空间 位 置 关 系 以 及 相对 速度 将 


它们 组 合成 最 终 模拟 使 用 的 IC XH 


准 。 如 D96 IF AHA 


EF, ABA, Xt M31 的 相关 测量 结果 就 显得 尤为 重要 。 

对 M31 的 距离 以 及 视 向 速度 ， 人 们 很 早 就 开始 研究 。 在 D96 中 ，Dubinski 等 人 选用 
了 700 kpc 的 星系 距离 以 及 -130lm s-! 的 径 向 速度 。 但 根据 M31 自行 计算 的 切 向 速度 直 
到 2012 年 才 被 正式 测定 ”"， 所 以 在 早期 的 工作 中 切 向 速度 的 选择 可 能 会 参照 某 些 特别 的 标 


E 粹 进行 银河 系 与 M31 并 合 模拟 的 工作 ， 主 要 目的 是 研究 潮汐 尾 与 暗 物 
质 野 间 的 关联 ， 所 以 通过 拟 合 潮汐 尾 的 形态 限制 并 设置 了 M31 的 切 向 速度 (20km . s- 和 


Ea 3 


该 


26 km - s~') 
置 了 相应 的 倾角 。 


2005 Æ, McConnachie 等 人 加 与 Ribas 等 人 四 研究 发 现 ， 氏 


作 已 经 考虑 到 星系 盘 平 面 与 轨道 平面 不 一 致 ， 给 银河 系 与 M31 设 


|5 M31 的 距离 大 约 


xy" 
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为 780kpc， 这 一 数值 被 C08 KA. Ait C08 的 开始 时 刻 并 不 是 现在 ， 
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时 刻 开 始 ， 当 时 银河 系 与 M31 的 距离 更 远 。 当 两 者 距离 演化 至 780 kpe, 


“初始 时 刻 ”， 而 根据 模拟 中 的 数据 ， 此 刻 M31 相对 银河 


向 速度 为 132 km . s-1， 这 也 是 对 于 M31 初始 速度 的 一 种 独特 设 定 方式 ， 


系 的 速度 和 倾角 并 不 完全 符合 实际 。 另 外 ，C08 的 最 终 IC 文件 


系 的 径 向 速度 为 -135km .s-1， 切 


其 模拟 是 从 更 早 的 
此 时 被 定义 为 新 的 


不 过 显然 会 导致 星 


还 包含 质量 M = 2.6 x 10 Mo M IGM. 


2012 年 ，M31 的 自行 终于 被 首次 测定 


中 除了 包含 银河 系 与 M31， 


" HST 得 到 的 相关 数据 "为 : M31 的 距离 


= (770 + 40) kpc， 相 对 银河 系 的 径 向 速度 Vind = (一 109.2 土 4.4)km.s !、 切 向 速度 


Vien = 17.0km:s-! (Æ 1o 的 置信 区 间 内 : Vien < 34.3km - s71). 


模拟 ， 并 且 在 模拟 中 加 入 了 M33， 是 一 个 三 星系 并 合 模 拟 。 
2019 年 ，van der Marel 等 人 "在 V19 中 利用 Gaia DR2 的 数据 给 出 M31 的 最 新 自 
行 结果 。 经 换算 ， 在 68% 的 置信 和 度 上 ，M31 切 向 速度 Vianpr2 = 1331km. st, 4A 


2012 aE HST 给 出 的 切 向 速度 Vian, HST = 36133 km . got, 
5723 km- sl。 工作 V19 最 后 便 是 利 


拟 还 考虑 了 M33 和 大 麦哲伦 云 ， 是 一 个 四 星系 并 合 模拟 。 


V12 采用 了 上 述 结果 进行 


加 权 平 均 的 结果 为 Vian pr24HsT = 
该 最 新 数据 进行 银河 系 与 M31 的 N 体 模拟 ， 此 模 


M31 最 新 切 向 速度 的 置信 下 限 和 置信 上 限 分 别 为 26km:s-! 和 92km.s-!， 置 信 区 间 


较 宽 ， 因 此 S20 固定 径 向 速度 (V. = 一 115.7 km s71) 


改变 切 向 速度 ， 并 进行 多 组 模拟 


图 四 展示 了 S20 中 银河 系 与 M31 的 初始 状态 。 然 而 ，S20 的 设置 略 有 不 足 ， 因为 人 类 观 


测 到 的 M31 视 向 速度 与 自行 是 在 太阳 系 内 观测 的 ， 而 进行 模拟 时 M31 
速度 是 相对 于 银河 系 (中 心 ) 的 ， 太 阳 系 的 位 置 与 银 心 的 位 置 存在 一 定 人 


的 径 向 速度 与 切 向 


局 差 ， 所 以 ,“ 自 行 ” 


对 应 的 方向 其 实 并 不 是 精确 的 “ 切 向 ”， 其 还 提供 了 极 小 部 分 M31 相对 于 银河 系 的 径 向 速 


度 ;“ 视 向 速度 ”对 应 的 方向 也 不 是 精确 的 “ 径 向 ”， 其 还 提供 了 一 小 部 分 M31 相对 于 银河 
能 完全 固定 径 向 速度 而 改变 切 向 速度 ， 同 时 也 不 能 仅仅 在 一 个 方向 
切 向 速度 矢量 应 当 不 是 在 一 条 直线 上 而 是 在 一 个 平面 内 变化 。 当 然 这 些 


系 的 切 向 速度 。 因 此 ， 不 


改变 切 向 速度 数值 


~ 


x 


“不 足 ”对 模拟 结果 影响 


的 影响 ， 只 是 从 更 严格 


注 : 此 图 为 S20 中 的 图 


2， 


y 


M31 


S 


Lilt, S20 的 这 一 操作 也 仅仅 是 为 了 研究 不 同 切 向 速度 对 并 合 时 间 
的 意义 上 来 讲 ， 这 些 设置 可 以 改进 


展示 了 模拟 的 初始 状态 下 银河 系 与 M31 的 位 置 、 速 度 矢 量 与 自转 角 动 量 矢量 。 


1 $20 中 银河 系 与 M31 的 初始 空间 关系 四 


改进 的 方法 可 以 参照 V12™”。 因 为 其 包含 3 个 星系 ， 显 然 无 法 简单 地 用 径 向 速度 和 切 
向 速度 将 银河 系 、M31、M33 进行 组 合 ， 必 须要 在 一 个 三 维 坐标 系 中 给 出 每 个 星系 具体 的 
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坐标 和 速度 矢量 ， 该 三 维 坐 标 系 便 


系 (Galactic coordinate system， 原 点 为 太阳 ),， 4 


是 银 心 参考 系 (Galactocentric frame)。 不 同 于 


41 卷 


坐标 


银 道 


心 参考 系 的 原点 为 银河 系 中 心 ， 因 


此 银 


河 系 的 初始 坐标 为 (0,0,0)kpc， 初 始 速度 为 (0,0,0) km : s-1。 银 心 参考 系 是 一 个 左手 系 ， 


其 中 zz 轴 指 向 太阳 ， 
度 、 赤 经 方向 


z 轴 指 向 北 银 极 。 我 们 可 以 将 
自行 (pmra)、 杰 纬 方向 自行 (pmdec) 转化 为 外 
速度 V,、V,、V.。 转 化 时 需要 根据 太阳 相对 于 银 心 的 位 置 以 及 速度 进行 修正 ， 


N 


hi A 


ERRIZ JRA EBR WER 


心 参考 系 下 的 坐标 x、y、z 与 
目前 普遍 


接受 的 数据 为， 太阳 与 银 心 的 距离 7 = 8.29kpc， 太 阳 所 在 半径 处 银河 系 的 旋转 速度 


Vo = 239km:s-!:， 太 阳 
PE 


出 每 个 星系 在 该 坐标 系 下 的 位 置 与 速度 ， 从 而 将 孤立 


4 星系 并 合 N 体 模拟 


最 终 的 IC 文件 构建 完成 后 ， 正 式 的 模拟 便 可 以 i 
体 模 拟 软 件 系 列 ， 同 时 附 上 近 些 年 来 太 


模拟 软件 。 本 章 将 主要 介绍 一 个 N 
作 中 针对 某 些 情形 积累 的 一 些 经 验 。 


4.1 GADGET 系列 


=I 


的 本 动 速度 (U,V,W)o = (11.1, 12.24, 7.25) km : s-1。 
合 模拟 工作 需要 选 定 一 个 合适 的 多 


4 


行 


总 而 言 之 ， 星 


E 标 系 ， 根 据 观测 者 位 置 与 该 坐标 系 原点 的 差异 ， 计 算 
系 的 IC 文件 正确 并 合 。 


— th A> 


款 合适 的 N 体 
究 人 员 在 模拟 工 


， 这 时 需要 选择 


NN 体 模拟 界 最 负 盛 名 的 当 属 Springe 等 人 开发 的 GADGET (GAlaxies with Dark 
matter and Gas intEracT)™ 系列 软件 ， 目 前 使 用 最 广泛 


的 GADGET-2™”， 下 面 对 其 进行 简要 介绍 。 


GADGET-2 是 一 个 能 够 在 具 


体 


模拟 的 核心 算法 为 TreePM 算法 


的 版 本 是 Springel 在 2005 年 发 布 


分 布 式 内 存 的 大 规模 3 
气体 动力 学 ) 的 开源 代码 ， 可 以 用 于 牛顿 力学 或 大 尺度 上 的 5 


行 计算 机 上 进 


宇宙 学 模拟 。 


行 N 体 模拟 (包括 
GADGET-2 中 N 


(也 用 到 了 SPH， 全 称 为 光滑 粒子 流体 动力 学 ， 其 与 气 


体 动力 学 的 计算 相关 ， 这 里 不 多 资 述 )，TreePM 算法 结合 了 Tree 与 PM (particle mesh, 


粒子 网 格 )， 不 同 于 PP 算法 ( 即 “粒子 -粒子 ”算法 ， 这 种 算法 是 最 基础 的 ， 直 接 求解 每 个 


N 


粒子 与 


他 所 有 粒子 间 的 相互 作用 ， 相 应 的 时 间 成 本 很 高 ， 与 粒子 数 的 平方 成 正比 )， 它 通 


过 构建 Tree 划分 粒子 ， 根 据 Tree 节点 判断 远近 将 粒子 间 的 作用 力 分 为 近 场 和 远 场 ， 近 场 


可 以 直接 计算 ， 远 场 则 采 
O(Nlg N) (N 为 粒子 数 )。 


在 利用 GADGET-2 进行 模拟 前 ， 使 月 


和 输出 文件 的 名 称 、 


JPM 算法 。 KA 


本 宇宙 学 参数 (对 于 牛顿 力学 的 模拟 3 
响 )， 计 算 机 相关 参数 (如 允许 运行 的 时 间 上 限 、 允 六 


过 相当 多 的 考量 ， 它 也 是 N 
说 明 。 


FE 大 大 提高 了 模拟 效率 ， 运 行 时 间 从 O(N?) 降 


F 


在 GADGET-2 的 使 用 者 手册 里 都 有 简单 介绍 与 参考 值 ， 但 其 中 软化 长 度 的 设 定 则 


T 


者 需要 在 参数 文件 中 填写 相应 的 参数 : IC 文件 
路 径 、 格 式 ， 模 拟 快照 和 数据 文本 的 输出 频率 ， 模 拟 的 时 间 跨 度 ， 基 
无 影响 ) 和 SPH 参数 (对 于 无 气体 粒子 的 模拟 并 无 影 
使 用 的 内 存 )， 软 化 长 度 等 。 所 有 参数 


i=} 
是 要 经 


本 模拟 中 极为 重要 的 一 个 参数 ， 本 文 的 4.2 et 


N 


设 定 作出 
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开始 模拟 后 ， 设 定 的 路 径 下 会 输出 一 系列 记录 模拟 过 程 的 文本 文件 (其 中 大 部 分 用 处 不 
大 ， 一 般 仅 在 模拟 出 问题 时 才 阅 读 这 些 文本 记录 以 进行 检查 ) 和 快照 文件 ， 编 写 代码 对 快照 
文件 进行 读 取 与 分 析 即 可 获得 数据 与 可 视 化 的 结果 。 


C08 选用 的 模拟 


学 元 


如 S20 中 使 月 


4.2 
4.2.1 


M31 


球 ) 粒子 的 


相等 


聚集 


A. AR 


Hi) HiGPUs™ 


软件 即 为 GADGET-2, 
GADGET 系列 的 最 新 版 软件 是 GADGET-4”， 它 在 代码 语言 、 
界面 等 方面 进行 了 优化 。GADGET 系列 软件 还 有 一 些 “ 衍 生 
GADGET-3 修改 而 来 的 GIZMO, GIZMO 的 功能 非常 之 广 ， 
射 、 黑 洞 吸 积 、 宇 


而 V12 选 月 


oo» 
HH » 


， 这 里 不 


N 体 模拟 常见 问题 处 理 


粒子 质量 的 设 定 


给 所 有 粒子 设置 相同 的 质量 


这 样 的 星系 ， 


MS, BORER 


在 中 心 ， 


4.2.2 星系 中 心 的 标定 


当 计 算 银 河 系 中 心 与 M31 4 


法 是 
黑洞 
的 质 
质量 


子 质量 ， 如 S20 中 就 将 银河 系 与 M31 的 中 心 黑洞 质量 设置 为 10? Mo 
然 偏 离 实 际 ， 但 在 并 合 发 生前 对 星系 尺度 的 动力 学 影响 甚 微 ， 如 原文 中 所 言 : 


距离 ， 不 过 严格 意义 上 来 讲 星系 3 
> 密度 最 高 点 或 者 势能 最 低 点 标定 中 心 ， 
Re AM aah TA 


月 


(如 C08) 似乎 是 一 
KAZET 90% 以 上 的 质量 ， 如 果 上 暗物质 晕 与 恒星 (包括 
质量 相等 ， 那 么 在 总 粒子 数 较 小 的 情况 下 ， 可 见 物 质粒 子 数 就 会 很 
究 者 对 于 某 些 结构 的 判断 和 认 知 。 所 以 ， 在 一 些 工 作 (如 V12) 中 不 同 种 类 粒子 的 质 
究 者 所 关注 的 部 分 则 会 
轻 粒 子 则 可 能 会 被 “弹出 ”， 
所 以 粒子 质量 也 不 宜 相差 过 大 


Ab E, 


He Bx 


量 可 
Æ H] 


a. 


4.2.3 


引力 


于 小 ; 


能 已 经 很 大 (比如 


这 时 便 不 能 将 


个 自然 的 选择 ， 


但 注意 到 对 


月 另 一 个 版 本 GADGET-3, 404 
软件 功能 


、 文 件 逻 辑 与 
比如 Hopkins 根据 


可 以 在 模拟 中 考虑 磁场、 化 
宙 线 等 。 除 了 GADGET 系列 以 外 ， 也 有 其 他 N 体 模拟 软件 ， 
再 介绍 。 


于 银河 系 与 


于 是 ， 


氏 点 。 但 由 于 计算 机 模拟 和 现实 情 
105 个 粒子 模拟 质量 为 1 x 10! Mo 
设置 为 实际 值 ， ny 


E. 
= 


黑洞 质量 


E 每 个 星系 中 心 置 入 大 质量 


星系 盘 与 核 
了 限 ， 影 响 研 
量 并 不 


由 更 多 粒子 组 成 。 不 过 在 引力 作用 下 ， 重 粒子 更 易 
系统 的 密度 分 布 、 稳 定性 与 演化 轨迹 将 受到 
， 有 具体 设置 视 所 研究 的 问题 而 定 。 


影响 ’ 


FP 心 的 距离 时 ， 可 以 计算 两 个 星系 的 质心 从 而 研究 质心 间 的 
不 是 对 称 的 结构 ， 质 心 并 不 等 同 于 中 心 ， 
可 以 大 


所 以 更 自然 的 方 
黑洞， 


存在 区 


区 别 ， 其 中 普通 粒子 


的 多 


河 系 时 ， 


每 个 粒子 的 


应 当 


软化 长 度 的 设置 


量 黑洞 只 是 星系 中 ， ÙD“ 被 动 的 客人 ”。 


软化 长 度 e (文中 采 ) 


发 散 的 数值 效应 ; 
Ay Fy MER 


计算 机 模拟 中 ， 
TY) 


星 的 


穿 过 彼此 而 不 是 碰撞 绚 


集合 ， 恒 星 碰 撞 截 


A AAT T FHE 


AJ 


的 单位 为 kpe) 是 N 体 模 拟 工作 中 的 重要 
当 两 个 粒子 距离 很 近 
HEN: N 体 模 拟 中 的 粒子 代表 的 不 是 
概率 极 低 ， Hi 
于 。 软 化 长 度 的 设置 有 一 个 经 验 公 式 : 


EAI, mi 


ag”, 


使 之 显著 大 于 普通 粒 


这 一 设置 虽 


时 ， pla 


— 


并 合 前 超大 质 


一 方面 避免 了 
中 的 分 母 值 过 


而 是 上 万 颗 恒 
星 交会 时 应 该 是 


(10) 
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KA, M 为 所 研究 的 单个 粒子 的 质量 。 经 分 析 ， 这 一 经 验 公 式 适 用 于 一 般 的 星系 并 合 模 
拟 ， 如 银河 系 与 M31 的 并 合 。 
软化 长 度 也 可 以 理解 为 模拟 的 空间 分 辩 率 ， 在 软化 长 度 尺度 下 的 动力 学 过 程 都 是 不 可 
靠 的 ， 存 在 非 物理 的 数值 效应 。 但 软化 长 度 并 不 是 越 小 越 好 ， 一 方面 ， GADGET-2 中 无 磁 
撞 粒 子 的 步 长 判 据 为 : 
2ne 
| "n 


其 中 ，At EK, Atmax 是 人 为 设 定 的 最 大 步 长 ，7 是 与 模拟 精度 相关 的 参数 ，a 是 粒子 
的 加 速度 。 因 此 ， 软 化 长 度 减 小 时 ， 模 拟 所 花 的 时 间 将 增加 。 另 一 方面 ， 如 前 文 所 提 ， 由 于 
粒子 数 有 限 ， 计 算 机 模拟 中 每 个 粒子 的 质量 其 实 很 大 ， 其 势 阱 也 非常 陡峭 ， 软 化 长 度 的 引入 
恰恰 使 引力 势 的 轮廓 更 加 平滑 ， 在 某 种 意义 上 能 较 好 地 符合 实际 情况 。 

软化 长 度 的 经 验 公式 也 不 唯一 。Rodionov 和 Sotnikova 给 出 依赖 于 粒子 数 的 经 验 
公式 : 


At = min 


os 


e 一 0.58NV -026 ， (12) 


其 中 ，N 为 模拟 中 的 总 粒子 数 ， 适 用 范围 为 30 < N < 300000. van den Bosch 和 Ogiya™ Jill 
针对 特定 的 暗物质 晤 粒子 给 出 经 验 公式 ; 


N.\ 73 
= 0.05 ( F2 ) , (13) 


RP, Ny 为 粒子 数 。 但 无 论 何 种 经 验 公式 都 只 是 一 个 参考 ， 最 终 的 设置 还 是 要 看 具体 的 模 
拟 需 求 ， 事 实 上 ， 许 多 工作 中 软化 长 度 的 设置 与 这 些 经 验 公式 相去 甚 远 。 


5 星系 并 合 模拟 结 


前 文 主要 选取 了 横 跨 25a 的 5 项 银河 系 与 仙女 座 星系 并 合 模拟 代表 性 工作 : D96、 
C08、V12、V19、S20。 现 将 这 些 工作 的 结果 进行 展示 。 

D96 并 不 是 专门 研究 银河 系 与 M31 并 合 的 工作 ， 所 以 给 出 的 结论 也 主要 着 眼 于 其 主 
题 一 一 星系 的 潮汐 尾 : 小 质量 的 银河 系 与 M31 模型 在 首次 交会 后 的 很 短 时 间 内 就 完成 并 
合 ， 并 合 时 刻 为 47Ga， 且 交会 过 程 中 会 产生 延伸 至 很 远 的 明显 潮汐 尾 。 大 质量 的 银河 系 与 
M31 模型 ， 由 于 加 速度 比 小 质量 模型 大 ， 第 一 次 交会 时 相对 速度 很 大 ， 只 产生 很 短 的 潮汐 
尾 并 且 潮 汐 尾 很 快 就 消失 ， 接 着 两 星系 分 离 至 反方 向 后 又 重新 接近 ;第 二 次 交会 时 也 未 能 
生 明显 的 潮汐 尾 ， 随 后 两 者 于 大 约 4.4 Ga 时 并 合 ” (模拟 的 可 视 化 结果 如 图 回 所 示 )。 

C08 的 模拟 结果 (星系 的 距离 -时 间 曲 线 见 图 团 显示 : 银河 系 与 M31 将 于 5Ga 以 内 完 
成 并 合 (此 时 太阳 尚未 死亡 )， 最 终 形 成 一 个 椭圆 星系 。 在 两 星系 第 一 次 交会 时 ， 太 阳 有 可 外 
被 抛射 到 一 个 潮汐 尾 中 ， 且 之 后 仍 有 较 大 概率 在 该 潮汐 尾 上 停留 ， 另外 太阳 也 有 较 小 的 可 外 
性 到 达 更 接近 M31 主体 的 区 域 ”。 
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iE: 此 图 为 D96 中 的 区 


13， 分 别 展示 了 小 质量 (上 ) 与 大 质量 (下 ) 模型 中 银河 系 与 M31 的 并 合 过 程 。 模 拟 中 


所 选择 的 切 向 速度 使 得 两 星系 产生 近 距 离 且 顺 行 的 交会 ( 顺 行 并 合 时 ， 盘 上 的 恒星 处 于 共振 态 ， 靠 近 扰动 体 一 侧 


的 恒星 被 持续 拉 进 ， 远离 扰动 体 的 恒星 被 持续 推出 一 )， 从 而 使 星系 盘 内 物质 受到 尽 可 能 大 的 作用 ， 最 终 小 


质量 模型 产生 了 较 明 显 的 潮汐 尾 ， 


大 质量 模型 则 未 能 成 功 产生 。 


2 D6 中 银河 系 与 M31 Het 


MW-M31 距 离 /Mpc 


1.4 鹃 T T T ™ T 


t/Ga 


YE: 此 图 为 C08 中 的 医 


5 


Ras SAR AG M31 距离 ( 纵 坐 标 ) 随时 间 ( 横 坐 标 ) 的 变化 曲线 。 水 平 虚线 对 应 的 


SERS 780kpc, PIE 
—4.7 Ga. 


Fan 


图 3 


星系 间 的 距离 ， 其 在 模拟 中 对 应 的 时 刻 约 为 4.7 Ga， 故 该 模拟 的 零 时 刻 相当 


Cos 中 银河 系 与 M31 的 距离 -时 间 图 像 
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V12 则 表明 ， 银 河 系 与 M31 将 于 3.87+032 Ga 时 第 一 次 交会 ， 近 心 距 为 31.0 
两 者 将 于 5.8615) Ga 时 完成 并 合 ( 见 图 团 。 数 值 存在 上 下 限 


自行 的 误差 范围 ， 采 用 了 “ 半 解 析 ” 的 蒙特 卡 洛 模拟 方法 。 对 了 
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十 38.0 
“i9.8 KPC 


[ 作 考 虑 了 M31 
F M31 和 M33， 两 者 第 一 次 


交会 的 时 间 为 0.857+0 呈 Ga， 近 心 距 为 80.8+43kpc。 对 于 银河 系 与 M33， 两 者 第 一 次 交会 


的 时 间 为 3.70+9.2 Ga， 近 心 距 为 176.07?399 kpc。 最 可 能 的 结果 是 银河 系 与 M31 并 合 后 ， 
M33 在 周围 轨道 上 绕 转 ， 稍 后 再 落 入 其 中 ; 但 也 有 M33 在 与 外 
落 入 的 可 能 ， 这 样 M33 将 先 于 M31 与 银河 系 并 合 ; Ab, M33 


次 交会 时 便 直接 
时 或 永久 地 被 


抛 出 本 星系 群 。 银 河 系 与 M31 的 并 合 产物 将 是 


地 方 ， 且 未 来 太阳 可 能 穿越 M33 。 


800 
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图 为 V12 中 的 图 3， 展 示 了 主 模拟 中 银河 系 、M31 5 M33 三 者 之 间距 离 (AA 


ik: 
变化 曲线 


4 V12 中 银河 系 、M31、M33 之 间 的 距离 -时 间 图 像 四 


V19 中 的 N 体 模 拟 并 不 关注 两 个 星系 最 终 的 并 合 ， 而 主要 着 
近 心 距 。 其 指出 ， 银 河 系 与 M31 BFA 4.5 Ga NRA 


第 一 次 交会 ， 近 心中 


而 考虑 M33 和 大 麦哲伦 云 的 影响 后 ， 首 次 交会 时 间 ] 


75 kpc ë 


心 更 远 的 


& 标 ) 的 


眼 于 首次 交会 的 时 间 与 
E 约 为 130 kpe; 
推迟 了 大 约 1Ga， 近 心 距 减 小 至 约 


S20 给 出 的 银河 系 与 M31 首次 交会 时 间 为 4.3Ga， 但 并 合 将 在 10 Ga 之 后 完成 ， 这 一 
时 间 比 前 面 的 工作 要 晚 很 多 ， 由 较 高 的 M31 切 向 速度 与 较 低 的 IGM 密度 共同 导致 ， 图 回 
展示 了 一 部 分 结果 。 该 工作 还 利用 ARWYV 代码 “研究 了 星系 并 合 后 两 个 超大 质量 黑洞 的 


行为 ， 指 出 两 个 黑洞 将 于 星系 并 合 后 的 16.6 Ma 内 并 合 。 另 外 该 2 


星系 的 暗物质 晕 半 径 、IGM 密度 、 模 拟 所 
洞 并 合 过 程 中 的 引力 波 辐射 。 


E: 
hu 


出 并 合 时 间 ， 所 以 此 处 未 展示 之 。 


JAF BOX HER 


HY M31 的 切 向 速 
并 计算 了 


以 上 工作 的 主要 参数 和 结果 见 表 区 由 于 V19 并 不 


整 的 并 合 模拟 工作 ， 未 给 
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YE: 此 图 为 S20 中 图 5 的 一 部 分 ， 展 示 了 切 向 速度 Vi = 57 km . s-! 的 情况 下 ， 不 同 IGM 密度 对 银河 系 与 
M31 并 合 时 间 的 影响 。 横 坐标 表示 时 间 ， 纵 坐标 表示 距离 。 


5 ”S20 中 银河 系 与 M31 的 距离 -时 间 图 像 四 


#2 D96, C08, V12 S20 的 主要 参数 与 结果 


参数 D96( 小 ) ”D96( 大 ) C08 V12 S20 
Mmw / 10° Mo 0.5 1-7 1.0 1.5 1.0 
Mms: / 10’? Mo 1.1 3.6 1.6 1.5 1.6 
Dinitial / kpc 700 700 780 770 780 
V; /km- s7! —130 —130 一 135 一 109.2 一 115.7 
V; /km-s7? 20 26 132 27.7 57 
IGM 含量 0 0 2.6 x 10'?Mo 0 10-26 kg - m~’? 
Tm / Ga 4.7 4.4 x5 5.86 >10 
HE: D96 的 小 质量 模型 和 大 质量 模型 无 主 次 之 分 ， 所 以 此 表 都 进行 了 展示 。 表 第 2 一 8 行 分 别 指 银河 


系 质量 、M31 质量 、 银 河 系 与 M31 初始 距离 、M31 径 向 速度 、M31 切 向 速度 、IGM 含量 (C08 
从 质量 角度 定义 、S20 从 密度 角度 定义 )、 并 合 完成 的 时 间 。 


O 


一 般 来 说 ， 银 河 系 与 M31 的 总 质量 越 大 、 初 始 距离 越 近 、 径 向 速度 越 大 ， 则 并 合 发 生 
的 时 间 越 早 ， 而 并 合 时 间 点 会 随 着 切 向 速度 的 增加 而 显著 推迟 ， 另 外 ， 当 IGM 的 数量 相当 
可 观 时 ， 其 提供 的 摩擦 会 使 得 星系 并 合 的 时 间 点 大 大 提前 ， 两 星系 于 首次 交会 后 不 久 即 发 生 


a 
H o 


6 总 结 与 展望 


星系 并 合 N 体 模 拟 工作 的 操作 流程 (如 图 四 所 示 ) 为 : (1) 构建 星系 模型 (确定 星系 质量 
等 参数 和 选择 星系 各 成 分 的 分 布 轮廓 )，(2) 生成 初始 条 件 (IC) 文件 (生成 孤立 星系 IC 文件 
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及 根据 相对 位 置 和 相对 速度 将 多 个 IC 文件 组 合成 一 个 新 IC 文件 ); (3) 正式 运行 模拟 ，(4) 


分 析 模 拟 结果 。 


构建 星系 模型 | 


emo | 


分 析 模 拟 结果 读 取 模拟 生成 的 snapshot， 进 行 数据 处 至 


6 ”星系 并 合 N 体 模拟 工作 操作 流程 


确定 基本 参数 : MB. HA. BRE. H 
选择 分 布 轮廓 : 决定 暗 晕 、 
生成 孤立 星系 IC 文件 

根据 各 星系 的 位 置 与 速度 将 IC 文件 组 合 


星系 盘 、 核 球 等 各 成 分 的 质量 分 布 


正式 运行 模拟 选择 合适 的 模拟 软件 ， 设 置 模拟 时 间 以 及 技术 性 参数 


与 可 视 化 


本 文 以 上 述 流程 为 线索 ， 介 绍 了 几 十 年 来 银河 系 与 M31 
选取 了 近 25 年 间 的 5 篇 代表 性 工作 (D96、C08、V12、V19、S20) 进行 比较 。 


开 合 模拟 工作 的 研究 进展 ， 并 


就 星系 质量 而 言 ， 早 期 人 们 的 估计 误差 范围 较 大 ， 如 今 普 裔 接受 的 结论 是 : 银河系 


的 质量 约 为 1.0 x 1012 Mo 或 者 1.5 x 10! Mo, M31 的 质量 约 为 1.5 x 1012 Mo 或 者 
2.0 x 10% Mo (但 也 有 最 新 结果 显示 其 可 能 接近 3.0 x 1012 Me)。 这 些 质量 在 星系 内 部 的 分 


布 则 可 以 用 一 些 轮廓 来 描述 ， 暗 物质 晕 常 


等 进行 描述 。 


常用 指数 轮廓 或 Miyamoto-Nagai 得 等 进行 描述 ， 核 球 常 


内 星系 的 空间 位 置 以 及 速度 关系 ， 可 以 将 
出 三 维 空间 中 具体 的 图 像 。 


GADGET 系列 软件 是 N 体 模 拟 界 最 著名 的 软件 ， 近 20 ER GADGET 发 行 了 不 少 新 
版 本 与 衍生 品 。 当 把 初始 条 件 文件 输入 N 体 模拟 软件 后 ， 将 会 有 一 系列 的 快照 文件 输出 ， 


读 取 与 分 析 这 些 快 照 文件 ， 我 们 将 一 览 星 系 


并 合 的 历程 。 


| Hernquist 轮廓 或 NEW 轮廓 进行 描述 ， 星 系 盘 
] Hernquist 轮廓 或 Plummer 球 


GALIC、DICE、MakeDisk 等 初始 条 件 生 成 器 的 出 现 ， 使 得 在 计算 机 中 创建 星系 粒子 
模型 的 操作 更 方便 ， 而 根据 星系 之 间 的 相对 位 置 和 相对 速度 将 各 个 星系 粒子 模型 进行 组 合 ， 
才能 够 得 到 N 体 模 拟 所 要 用 的 最 终 初始 条 件 文件 。 近 几 十 年 来 ， 人 们 对 M31 的 位 置 和 径 向 
速度 认 知 的 改变 并 不 大 ， 如 今 测 得 M31 的 距离 约 为 770kpc， 径 向 速度 约 为 -109km . s-1。 
而 M31 的 自行 直到 2012 年 才 被 HST 首次 正式 测定 ， 又 于 2019 年 被 Gaia 再 次 测定 ， 将 两 
个 结果 加 权 平 均 并 换算 得 出 的 M31 切 向 速度 为 57+33 km . s-1。 为 了 更 准确 地 描述 本 星系 群 
星系 的 位 置 和 速度 转换 到 银 心 参考 系 下 ， 从 而 给 


本 文选 取 的 5 篇 代表 性 工作 ， 尽 管 初始 条 件 各 异 ， 但 银河 系 与 M31 的 首次 交会 时 间 均 
约 为 4Ga; 而 最 终 并 合 时间 却 因 模拟 参数 的 不 同 而 大 相 径 庭 。 这 些 工 作 也 并 非 只 研究 了 安 


观 上 的 并 合 结果 ， 还 关注 太阳 的 未 来 命运 、 
银河 系 与 M31 的 并 合 模拟 工作 已 经 取 


星系 中 心 黑洞 的 并 合 等 。 


得 了 相当 大 的 成 功 ， 但 依然 有 可 以 改进 的 地 方 : 


(1) 随 着 观测 能 力 的 进一步 提高 ，M31 的 自行 和 本 星系 群 的 质量 估计 有 了 新 的 结果 ， 这 需要 
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在 模拟 中 做 出 相应 的 改变 ，(2) 部 分 参数 的 误差 范围 仍然 较 大 ， 需 要 进行 覆盖 全 部 误差 范围 
的 模拟 工作 ; (3) 本 星系 群 中 除了 银河 系 与 M31 之 外 其 他 成 员 的 影响 也 不 可 忽略 ， 多 星系 模 
拟 可 能 成 为 未 来 的 一 种 趋势 ，(4) 包含 复杂 物理 过 程 的 模拟 越 来 越 热 门 ， 银 河 系 与 M31 未 来 
并 合 过 程 中 的 恒星 形成 、 黑 洞 生长 等 将 是 很 有 价值 的 课题 。 
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N-body Simulations 
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Abstract: The Local Group contains nearly one hundred galaxies of various sizes and 
shapes, including the Milky Way (MW), the Andromeda Galaxy (M31), the Triangulum 
Galaxy (M33), the Magellanic Clouds and so on. M31 is the most massive and brightest 
galaxy in the Local Group. Studies of the Local Group have significant implications for 
understanding the interactions and evolution of galaxies, as well as cosmology. According to 
the hierarchical clustering theory of the universe, galaxies grow through mergers. In addition 
to numerous satellite galaxy mergers in the Local Group, there may also be a future merger 
between the two large galaxies, the Milky Way and M31. Observations indeed indicate 
that M31 and the Milky Way are approaching each other and will merge in the future. 
As a powerful tool for studying galaxy mergers, N-body simulation has rapidly matured 
during the past half century and played an important role in our quantitative understanding 
of galaxy mergers. Thanks to the development of galaxy structure theories, observational 
capabilities and N-body simulation software, information about the Milky Way and M31 has 
become increasingly clear, thus our understanding of the merger between the Milky Way and 
M31 has become more detailed. This article presents a brief overview of N-body simulations 
of the merger between M31 and MW, with an emphasis on the simulation procedure, as well 


as a comparison of recent results. Valuable future directions are also discussed. 
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